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An  elastodynaraic  conservation  integral,  the  integral,  is  employed  to 
derive  boundary  integral  equations  for  crack  configurations,  in  a  direct  and 
natural  way.  These  equations  immediately  have  lower  order  singularities 
than  the  ones  obtained  in  the  conventional  manner  by  the  use  of  the  Betti - 
Rayleigh  reciprocity  relation.  This  is  an  important  advantage  for  the 
development  of  numerical  procedures  for  solving  the  BIE's,  and  for  an 
accurate  calculation  of  the  strains  and  stresses  at  internal  points  close  to 
the  crack  faces.  For  curved  cracks  of  arbitrary  shape  the  BIE's  presented 
here  have  simple  forms,  and  they  do  not  require  integi-ation  by  parts,  as  in 
the  conventional  formulation.  For  the  dynamic  case,  the  unknown  quantities 
are  the  crack  opening  displacements  and  their  derivatives  (dislocation 
densities),  while  for  the  static  case  only  the -dislocation  densities  appear 
in  the  formulation.  For  plane  cracks  the  boundary  integral  equations  reduce 
to  the  ones  obtained  by  other  investigators. 


i  1)1  :'f  ribut  I  on/ 

Avnilablllty  Codes 
r  lAvail  an4/or 


1 .  Introduction 

Boundary  integral  equations,  in  conjunction  with  the  boundary  element 
method,  provide  an  effective  numerical  technique  for  the  solution  of 
boundary  value  problems  in  solid  mechanics.  The  boundary  integral  equation 
method  (BIEM)  has  been  successfully  applied  to  a  wide  range  of  problems  in 
linear  and  nonlinear  elasticity.  Recent  developments  of  the  boundary 
integral  equation  method  have  been  concerned  with  applications  to 
elastostatic  and  elastodynamic  crack  analysis.  The  method  is  attractive  for 
crack  analysis,  because  the  semi-analytical  nature  of  the  BIEM  makes  it  easy 
to  take  into  account  the  singularities  at  the  crack  tips. 

The  conventional  BIE  formulation,  due  to  Rizzo  (1967)  and  Cruse  (1969) , 
is  based  on  the  Betti -Rayleigh  reciprocity  theorem  for  two  independent 
elastostatic  or  elastodynamic  states.  By  choosing  one  of  the  states  as  the 
unknown  field  and  the  other  as  the  basic  singular  solution  (the  Green's 
function) ,  a  representation  integral  for  the  displacement  components  can  be 
derived.  The  integral,  which  is  over  the  surface  of  the  crack  contains  the 
crack  opening  displacements  (the  displacement  jumps  across  the  crack  faces) 
and  derivatives  of  the  Green's  function  in  its  integrand.  Unfortunately,  a 
direct  limiting  process  on  the  representation  integral  for  the  displacements 
as  the  observation  point  approaches  a  crack  face,  gives  rise  to  a  degenerate 
set  of  BIE's,  as  shown  by  Cruse  (1975).  This  has  motivated  the  use  of 
representation  integrals  for  the  tractions,  and  their  corresponding  boundary 
integral  equations,  rather  than  displacement  BIE's.  Such  traction  BIE's 
are,  however,  highly  singular,  and  they  cannot  be  solved  directly  by 
numerical  methods.  To  circumvent  these  difficulties  several  approaches  have 
been  proposed,  see  for  example  the  papers  by  Cruse  (1975,1987),  Weaver 
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(1977),  Budiansky  and  Rice  (1979),  Schmerr  (1982),  Sladek  and  Sladek  (1984), 
Nishimura  and  Kobayashi  (1987),  Zhang  and  Achenbach  (1988)  and  Budreck  and 
Achenbach  (1988) .  Most  of  these  studies  first  reduce  the  higher  order 
singularities  to  integrable  ones,  and  then  solve  the  modified  BIE's 
numerically.  The  reduction  is  achieved  by  the  use  of  partial  integration. 
The  required  manipulations  are  reasonably  easy  for  simple  configurations 
such  as  3-D  planar  or  2-D  straight  cracks,  but  they  become  cumbersome  for 
curved  cracks.  Furthermore,  different  forms  of  the  regularized  BIE's  are 
obtained  through  the  non-unique  integration-by-parts  process,  though  they 
are  equivalent  (see  Cruse,  1987). 

In  this  paper  we  present  a  new  BIE  formulation  for  crack  analysis.  The 
motivation  for  this  study  is  the  paper  by  Hu  (1987) ,  who  proposed  a  novel 
way  to  obtain  a  new  type  of  BIE,  to  solve  elastostatic  boundary  value 
problems.  Hu's  formulation  is  based  on  the  conservation  integral  Jj^.  In 

the  present  paper  it  is  shown  that  Hu's  BIE's  are  especially  suited  for 
solving  crack  problems.  For  elastodynamic  problems,  the  integral  of 

elastostatics  is  generalized  to  time-harmonic  elastodynamics ,  and  the  result 

is  denoted  by  Boundary  integral  equations  are  then  derived  from  in  a 

direct  and  natural  manner,  for  arbitrary  crack  configurations.  The  BIE's 
that  are  obtained  are  immediately  less  singular  than  the  ones  of  the 
conventional  formulation,  and  they  do  not  require  additional  manipulation  in 
developing  numerical  solution  procedures.  The  BIE's  presented  in  this  paper 
nave  lel^tively  simple  forms,  and  they  reduce  to  those  for  elastostatics  by 
letting  -*  0  and  by  using  the  appropriate  stTtlc  Kelvin  solutions.  For 
planar  and  straight  cracks  the  results  agree  with  those  obtained  by  other 
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authors.  New  BIE's  derived  from  the  complementary  conservation  integral 

are  also  given,  but  these  equations  do  not  offer  advantages  over  the 
conventional  formulation. 

The  significance  of  (or  Jj^)  integral  as  a  relevant  crack- tip 

parameter  in  linear  and  nonlinear  fracture  mechanics  has  been  well 
established  (see  Moran  and  Shih,  1987).  The  present  paper  presents  a  novel 
application  of  such  path  independent  or  conservation  integrals  in 
elastodynamic  and  elastostatic  crack  analysis. 

2.  Problem  Statement  and  Conventional  BIE_Formulation 


A  crack  is  a  surface  of  displacement  discontinuity  when  external  loads 
are  applied  to  the  body.  The  faces  of  a  mathematical  crack  are 
infinitesimally  close  prior  to  loading,  and  they  do  not  interact  when  loads 
are  applied.  This  is  an  acceptable  approximation  for  real  cracks  whose 
faces  are  initially  sufficiently  separated  so  that  the  faces  will  not  touch 
when  the  body  is  disturbed. 

In  this  paper  we  consider  a  three-dimensional  (curved)  crack  of 
arbitrary  shape  which  is  contained  in  an  unbounded,  homogeneous,  isotropic, 
linearly  elastic  solid.  The  geometry  is  shown  in  Fig.  1.  The  solid  is 
subjected  to  time -harmonic  motion,  but  the  term  exp(-iwt)  has  been 
suppressed  throughout  the  analysis. 

The  stress  equations  of  motion  are  given  by  (see  Achenbach,  1973) 


o . .  .  +  pw^u.  —  0  , 

i-J.J  ^ 
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(2.1) 


where  defines  the  stress  components,  denotes  the  displacement 

components,  p  is  the  mass  density,  and  a>  is  the  angular  frequency.  In 
Eq.(2.1)  body  forces  are  not  considered  and  the  summation  convention  is 
implied.  In  the  linear  theory  the  strain  components  are  defined  as 


€ .  .  -  r  (u.  .  +  u.  . ) 
i-J  2  i,j  j,i' 


The  stress  and  strain  components  are  related  by  Hooke's  law 


(2.2) 


(2.3) 


where  are  elastic  constants  which  for  isotropic  materials  can  be 

written  as 


^ijki  ■  ^-^ij^ki  ^ii-jk^ 


(2.4) 


Here  A  and  u  are  Lame's  elastic  constants  and  5..  is  the  Kronecker  delta. 
The  tractions  vanish  on  the  faces  of  the  crack,  i.e.. 


f.-c7..n.-0  ,  xeA  , 

1-  i-J  J  - 


(2.5) 


where  A  =  +  A  .  For  a  scattering  problem  a"*"  is  the  insonified  side  of 

the  crack  and  A  is  the  shadow  side.  Also,  n^  is  the  unit  normal  vector  of 


The  total  fields  generated  by  the  interaction  of  an  incident  wave  with 
the  crack  can  be  written  as 


in  ,  sc 
u.  -  u.  +  u. 
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in  sc 
a .  .  —  a .  .  +  a .  . 

i-J 
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where  and  represent  the  incident  field  in  the  absence  of  the  crack, 
SC  sc 

and  define  the  scattered  field.  Both  the  total  fields  and  the 

partial  fields  satisfy  Eqs . (2 . 1) - (2 . 3) .  For  a  given  incident  field  the 
scattered  field  has  to  satisfy  the  boundary  conditions  on  the  faces  of  the 
crack,  Eq. (2 . 5) . 

Following  the  procedure  proposed  by  Rizzo-  (1967)  and  Cruse  (1969),  a 
representation  integral  for  the  scattered  displacement  can  be  obtained  by 
using  the  Betti-Rayleigh  reciprocal  theorem  and  the  fundamental  solution  due 
to  a  unit  time-harmonic  point  force.  For  a  3-dimensional  crack,  the 
representation  integral  can  be  written  as 


u^^(x)  =  <7^j^(x-y)Au^(y)njdA(y) ,  x  i  a"*"  . 


(2.7) 


Here  x  is  the  position  vector  of  the  observation  point,  aT.,  (x-v)  is  the 

ij  K  -  — 

stress  Green's  function  (Appendix  A),  and  AUj^(y)  defines  the  displacement 

jumps  (crack  opening  displacements)  across  the  faces  of  the  crack. 

As  shown  by  Cruse  (1975)  for  the  static  case,  Eq.(2.7)  will  lead  to  a 

degenerate  BIE  formulation  as  x  -*■  A^.  A  natural  remedy  for  this  difficulty 


is  to  use  the  representation  integral  for  the  traction  components,  which  can 
be  obtained  by  substituting  Eq.(2.7)  into  Hooke's  law  and  by  using 


-  cr^^n  .  The  result  is 

p  pq  q 
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Boundary  integral  equations  can  be  derived  from  Eq.(2.8)  by  letting  x  -» 

and  by  applying  the  boundary  conditions  (2.5).  The  system  of  boundary 
integral  equations  that  is  obtained  in  this  manner  is,  however, 
hypersingular  when  the  observation  point  x  and  the  source  point  ^  coincide, 

Q 

since  the  terms  behave  as  (Appendix  A) 
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(2.9) 


where  r  -  |x-;^| .  These  higher  order  singularities  prevent  a  reliable  direct 
numerical  solution  of  Eq.(2.8). 

To  overcome  these  difficulties  Budiansky  and  Rice  (1979)  used  partial 
integration  to  reduce  the  higher  order  singularities,  and  to  derive  a  system 
of  BIE's  for  a  flat  crack  in  the  plane  x^  -  0  (see  Fig.  1).  Regularization 
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procedures  have  also  been  proposed  by  Sladek  and  Sladek  (1984),  by  Nishimura 
and  Kobayashi  (1987),  and  by  Budreck  and  Achenbach  (1988). 


For  a  two-dimensional  crack  configuration,  analogous  formulations  have 
been  proposed  by  Tan  (1975),  by  Schmerr  (1982),  and  by  Zhang  and  Achenbach 
(1988).  The  corresponding  elastostatic  crack  analysis  using  BIE  methods  has 


been  presented  by  Cruse  (1975)  and  Weaver  (1977) .  A  comprehensive 
discussion  and  an  extensive  list  of  references  has  been  given  by  Cruse 
(1987). 

All  the  studies  mentioned  above  have  used  partial  integration  to  reduce 
the  higher  order  singularities  (2.9).  This  procedure  is  easily  implemented 
for  flat  or  straight  cracks,  but  it  becomes  quite  cumbersome  for  curved 
cracks.  In  this  paper,  we  will  present  a  new  BIE  formulation  which  follows 
very  naturally  from  a  path  independent  integral ,  and  which  has  lower  order 
singularities  than  the  ones  obtained  in  the  conventional  BIE  formulation. 

3.  The  Integral  and  Related  BIE's 
In  elastostatics  the  integral  has  the  form  (Eshelby,  1951,  Rice, 

1968) 


a .  .  u .  ,  )  n .  dS 


(3.1) 


where  S  is  the  surface  of  a  body  with  volume  V,  n^  is  the  outward  normal 


vector,  and  W  is  the  elastic  strain  energy  density 


W 


1 

r  a .  . 

2  i-J  ij 


(3.2) 


The  integral  which  vanishes  if  there  are  no  body  forces  and 

singularities  present  in  V,  is  usually  referred  to  as  a  path  independent 
integral  or  a  conservation  law.  The  application  of  the  component  as  a 

relevant  crack- tip  parameter  in  linear  and  nonlinear  fracture  mechanics  has 
been  well  established  (see  Moran  and  Shih,  1987).  It  can  easily  be  shown 
that  =•  0 ,  by  applying  the  divergence  theorem,  by  using 


aw _ 

‘^ij  ~  3i.  . 


(3.3) 


and  by  employing  the  equilibrium  equations  ^7^^  ^  -  0.  The  generalization  of 
to  time -harmonic  elastodynamics ,  which  is  denoted  by  can  be  written 


(3.4) 


where  L  is  the  kinetic  energy  density 


L  -  i  pii.u. 


1  2 
—  pcj-'u.u. 
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(3.5) 


Here  also  »  0,  under  the  same  assumptions  as  for  The  proof  is  again 

very  simple  if  we  apply  the  divergence  theorem,  use  Eq.(3.3)  and  employ  the 
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equations  of  motion  (2.1).  We  note  that  -  0  holds  for  any  time -harmonic 

elastodynamic  state  which  satisfies  Eqs . (2 . 1) - (2 . 3) . 

Now  let  us  consider  two  independent  time -harmonic  elastodynamic  states 
for  the  same  body: 


l-J  3-J 


(3.6) 


.  ,(2) 

i-J 


(3.7) 


These  states  satisfy  the  equations  of  motion  (2.1),  the  strain-displacement 
equation  (2.2)  and  Hooke's  law  (2.3).  By  virtue  of  linear  superposition, 
the  sum  of  (3.6)  and  (3.7) 


.  u<2)  .  +  e(2)  , 


U.  “  u.  +  u. 

1  r  I 


l-J  l-J  3-3 


a.  .  -  ,  (3.8) 
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also  satisfies  Eqs . ( 2 . 1) - ( 2 . 3) .  Substitution  of  (3.8)  into  (3.4)  yields 


J,  [u. ] 
k  1 


J  [u^^^ ]  +  J  [u^^^ 
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Clearly,  the  terms  1  and  ]  must  vanish  because  u^  and  u^ 

are  two  independent  elastodynamic  states.  Since  [u^]  -  0  we  therefore 
conclude  that 


;  .  a^!^^"^n.dS  -  0  .  (3.10) 

^^  m,nmn  i  i  jk  i,kLj  ijL.k^j 
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Equation  (3.10)  is  an  extension  of  the  two-state  conservation  integrals 
proposed  by  Chen  and  Shield  (1977)  for  elastostatics  (co  -»  0)  .  The  first 
state  is  now  taken  to  be  the  unknown  field 


'  L  '  LJ 


,  SC  SC, 

{u.  ,a. .] 
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(3.11) 


while  the  second  state  is  selected  as  the  fundamental  solution  due  to  a  unit 
point  force 


(u^2),,(2) 

1  LJ 


(3.12) 
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where  u. „  and  a...  are  3-D  time-harmonic  elastodynamic  Green's  functions 

li  Lji  ^ 

(see  Appendix  A) ,  and  a^  indicates  the  directions  of  the  applied  point 
force.  Application  of  (3.10)  to  the  surfaces  S,  and  (Fig. 2),  and  use 
of  Che  Eqs.(3.11)  and  (3.12)  results  in 
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I  i^i,(x;v)dS(v)  +  J  =  ^  I  -  0  • 
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where  S  is  the  surface  of  the  scatterer,  is  the  surface  of  a  sphere  of 

radius  5,  centered  at  x,  and  is  the  surface  of  a  sphere  with  radius  R, 

centered  at  x,  as  shown  in  Fig.  2.  The  surface  S  is  assumed  to  be  closed, 

regular  and  smooth.  The  small  sphere  (radius  S)  is  selected  to  exclude  the 

singularities  in  the  Green's  functions,  and  the  sphere  with  radius  R  must  be 

sufficiently  large  so  that  the  scatterer  S  and  the  sphere  S ,  are  included  in 

0 

it.  The  integrand  1 in  (3.13)  is  given  by 


SC 


(3.14) 


in  which  .x  represents  the  position  vector  of  the  observation  point,  and  ^ 

represents  the  position  vector  of  the  source  point,  respectively. 

After  elementary  calculations  the  second  integral  in  Eq.(3.13)  can  be 
evaluated  as 
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(3.15) 
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By  using  the  asymptotic  expansions  of  the  Green's  functions  for  large  ly| 

(Appendix  B)  the  last  integral  in  Eq.(3.13)  can  be  rewritten  in  the 
following  form 


A+2/1 


exp(ik^R) 
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Applying  the  Cauchy- Schwartz  inequality 


i;  f(ir)dS(jr)l  <  (/  dS(y))''(;  IfCyll^dSCy))"" 
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(3.1 


and  considering  the  following  elastodynamic  radiation  conditions  (see  Tan, 


1975,  Achenbach,  1982) 
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we  obtain 


Js  I^^(x:ir)dS(^)  -  0 
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as  R  -►  ® 


(3.20) 


Thus,  Eq.(3.13)  is  reduced  to 


sc  ,  . 


X  ^  S . 


(3.21) 


Substitution  of  Eq.(3,21)  into  Hooke's  law  leads  to  the  following 
representation  integral  for  the  traction  components  at  x 


-SC  ,  , 

fp  ^2) 


sc ,  ,  G 


sc ,  ^  G 


(3.22) 


Equation  (3.22)  is  a  generalization  of  Hu's  results  for  elastostatics 
(1987) . 


Application  of  Eq.(3.22)  to  a  3-D  crack  yields 


irSC  ,  , 

fp  (2)  - 


-  Au^  .  X  'I  A'^.  (3.23) 

where  Au.  are  the  crack  opening  displacements  and  Au^  are  their 
derivatives  with  respect  to  y^.  The  last  term  of  Eq.(3.22)  disappears 
because  of  the  continuity  of  across  Che  crack  faces.  BIE's  are 

obtained  by  letting  x  -»  A^  as 


%  (5) 


G  + 

A'^i  .  X  e  A  . 


(3.24) 


The  integral  of  (3.24)  is  understood  in  the  sense  of  Cauchy  principal 
values.  No  extra  discontinuity  terms  enter  (3.24)  as  x  -►  A^.  Equation 
(3.24)  is  valid  for  a  3-D  crack  of  arbitrary  shape.  The  corresponding  BIE's 


for  a  2-D  crack  in  plane  strain  and  anti-plane  strain  can  be  derived 
directly  from  (3.24).  For  plane  strain  we  obtain 


while  for  anti -plane  strain  we  find 

f^^Cx)  -  M  ^  1 

-  AUj  '  X  e  r'*'.  (3.26) 

Here  denotes  the  insonified  side  of  the  2-D  crack  (see  Fig.  3)  and  the 
superscript  "g"  represents  the  2-D  Green's  functions  (Appendix  A).  All 
integrals  of  (3.25)  and  (3.26)  are  understood  as  Cauchy  principal  values. 

The  BIE's  (3.24)  (as  well  as  (3.25)  and  (3.26))  have  the  advantage  over 
the  conventional  BIE  formulation,  i.e.,  Eq.(2.8),  that  no  higher  order 
singularities  appear.  The  unknown  boundary  quantities  in  the  new 
formulation  are  the  crack  opening  displacements  and  their  derivatives,  where 
the  latter  have  the  physical  meaning  of  dislocation  densities.  We  note  also 
that  the  procedure  in  deriving  (3.24)  is  very  natural,  and  no  elaborate 
manipulations,  such  as  integration  by  parts,  have  been  used.  The  BIE's 
stated  here  apply  also  to  elastostatic  crack  analysis  by  letting  u  -*  0 ,  and 
by  using  the  corresponding  elastostatic  fundamental  solution  (Kelvin 
solution).  For  the  static  case,  the  term  containing  Au^  disappears  in  the 
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BIE's,  and  the  only  unknown  quantities  are  the  dislocation  densities.  The 


new  formulation  allows  an  immediate  numerical  implementation.  When  Au^  has 

been  computed  equation  (2.7)  can  be  employed  to  calculate  the  displacement 
field. 

As  pointed  out  by  Carlsson  (1974),  Bui  (1977),  and  Moran  and  Shih 
(1987),  there  exists  a  complementary  integral  to  which  is  also  path 

independent.  For  time-harmonic  elastodynamics ,  the  complementary  integral 
may  be  stated  as 


3'.;  ((u-= 


dS  -  0  , 


(3.27) 


c  c 

where  W  and  L  are  defined  by  the  Legendre  transformation 


W^  -  <7.  .  -  W  , 

ij 


L  —  pii.u.  -  L  , 
11 


(3.28) 


(3.29) 


and  W  and  L  are  given  by  (3.2)  and  (3.5),  The  assumptions  in  deriving 
(3.27)  are  the  same  as  for  J^.  The  proof  of  Eq.(3.27)  can  be  performed 
directly  by  using  the  divergence  theorem,  by  considering  the  relation 


_  aw 

'  ij  "  d<7.  . 

i-J 


(3.30) 


and  by  employing  Eq.(2.1). 


Following  the  same  procedure  as  in  deriving  Eq.(3.22),  a  novel 
representation  integral  for  the  scattered  traction  components  is 

obtained  from  (3.27)  as 


-  c 


■  ^ij  .  5  ^  S.  (3.31) 

Corresponding  BIE's  can  be  derived  by  applying  Eq.(3.31)  to  the  crack  faces 
and  by  letting  x  -►  A,  taking  into  account  of  the  boundary  conditions  (2.5). 

Such  BIE's  are,  however,  again  highly  singular  due  to  the  presence  of  the 

Q 

term  •  Hence,  Eq.(3.31)  has  no  advantages  over  the  conventional 

formulation  given  by  Eq.(2.8). 

4 .  Examples 

In  this  section  we  will  apply  the  BIE's  (3.24)  (as  well  (3.25)  and 
(3.26)),  which  are  valid  for  arbitrary  shaped  cracks,  to  some  simple  cases. 
We  first  consider  a  flat  3-D  crack  in  an  unbounded  body  subjected  to  an 

incident  time -harmonic  wave.  The  crack  is  located  in  the  plane  x^  “  0“. 

Hence  n^  -  n2  -  0  and  n^  =  1.  The  BIE's  (3.24)  separate  into  two  decoupled 

equations : 

^723(x^,X2,0)  -  {  { (A+2,u)a^33(x-y)  -  Ao33^(x-y)  jAu^  ^(y)  - 

A  ’ 


-  (.\+2 pi) Ayi^(l)  )dAi2)  ,  5  «  A 


(4.1) 


in 

^3 


(x^,X2.0)  -  - 


(4.2) 


where  *^^3  stress  components  corresponding  to  the  incident  wave. 

We  note  that  Eq.(4.1)  is  for  the  normal  crack  opening  displacement  Au^  while 

Eq.(4.2)  is  for  the  transverse  crack  opening  displacements  Equations 

(4.1)  and  (4.2)  have  exactly  the  same  forms  as  those  derived  by  Budiansky 
and  Rice  (1979)  who  used  the  conventional  formulation  in  conjunction  with 
partial  integration. 

BIE's  for  a  flat  3-D  crack  under  static  surface  loading  ff22(x^,X2,0) 
and  a^^(x^,X2,0)  can  be  obtained  from  (4.1)  and  (4.2)  by  letting  w  -►  0  and 
by  employing  the  corresponding  static  fundamental  solutions .  The  result  is 


a33(x^,X2.0)  =  S  ^^3,0^2)  dA(y).  x  .  A+  , 


<7^3(Xi,X2.0)  -  8^.5!-i/)  ■^+  r2  ^  ^ 


(4.3) 


(4.4) 


Here  r  —  |x-^l ,  and  i/  denotes  Poisson's  ratio.  Equations  (4.3)  and  (4.4) 

are  identical  to  the  equations  stated  by  Weaver  (1977). 

Next,  we  consider  a  straight  2-D  crack  for  states  of  deformation  in 

+ 

plane  strain  and  anti-plane  strain.  The  crack  is  defined  by  X2  -  0  , 

|x^|  <  a.  For  the  case  of  plane  strain,  the  BIE's  (3.25)  become 


-pw2u?^(x-y)Au, (y, ) )dy^,  jx-l  <a  , 


O22(Xx.0)  -  J  ([(>  +  '  ^'^221  ^  ^'^2 , 1  ^^1^  ‘ 


■pu'^iX  +  2pt)u^^(x-i)^^2^yi))dy^,  |xj^l  <  a  , 


while  for  anti-plane  strain  Eq.(3.26)  takes  the  following  form 


(722(^1-0)  "4/  I  l^^l^  * 


|x^|  <  a  . 


(4.5) 


(4.6) 


(4.7) 


mm 


The  BIE's  (4.5)  and  (4.6)  have  been  derived  by  Tan  (1975)  via  Che 


conventional  formulation.  For  plane  strain  deformation,  the  BIE  for  the 


normal  crack  opening  displacement,  AU2  (Mode  I),  decouples  from  the  one  for 


the  transverse  crack  opening  displacement,  Au^^  (Mode  II).  For  a  2-D  crack 


under  static  loading  obtain  from 


(4.5)  and  (4.6) 


a  Au, 


cr^2(Xi.O)  -  2ir(l-i/)  x^-y^  ^^^1  ’ 


l^il  <  a 


(4.8) 


a  Au. 


ct22(Xi.O)  -  2^(1-;,)  ^I’Yi  ’  l^ll  -  ^  ’ 


(4.9) 


for  plane  strain,  and  from  (4.7) 


a  Au- 


a32(Xi,0)  -  -  dy^  . 


1x3^1  <  a  , 


(4.10 


for  anti-plane  strain.  Equations  (4.8)  -  (4.10)  are  integral  equations  for 


dislocation  densities,  and  they  are  again  well  known  (see  Mura,  1987). 


The  BIE's  presented  here  must,  in  general,  be  solved  numerically. 


Special  care  must  be  taken  in  the  numerical  implementation  to  account  for 


the  local  behavior  of  Au.  and  Au.  .  near  crack  edges,  and  for  the 

^  1-  •  J 


singularities  of  the  Green's  functions  at  x  -  For  2-D  cracks  subjected 


by  static  loading  the  method  developed  by  Erdogan  et  al.  (1973)  has  been 


frequently  used,  while  Zhang  and  Achenbach  (1988)  solved  the  modified  BIE's 


of  (4.5)  and  (4.6)  numerically  for  time-harmonic  wave  scattering  problems. 


For  a  flat  3-D  crack  numerical  methods  have  been  proposed  by  Polch  et  al . 


(1987)  for  the  static  case,  and  by  Nishimura  and  Kobayashi  (1987)  for  the 


dynamic  case. 


5 .  Concluding  Comments 


A  novel  application  of  an  elastodynamic  conservation  integral,  the  J, 


integral,  to  elastodynamic  and  elastostatic  crack  analysis  has  been 


presented.  Boundary  integral  equations  follow  from  in  a  direct  and 


natural  way.  These  equations  immediately  have  lower  order  singularities 


than  the  ones  obtained  in  the  conventional  manner  by  the  use  of  the  Betti - 


Rayleigh  reciprocity  integral.  This  is  an  important  advantage  for  the 


development  of  a  numerical  procedure  for  solving  the  BIE's,  and  for  an 


accurate  calculation  of  the  strains  and  stresses  at  internal  points  close  to 


the  crack  faces.  For  3-D  or  2-D  cracks  of  arbitrary  shapes  the  BIE's 


presented  here  have  simple  forms,  and  they  do  not  require  integration  by 


parts,  as  in  the  conventional  formulation.  In  the  dynamic  case,  the  unknown 


quantities  are  the  crack  opening  displacements  and  their  derivatives 


(dislocation  densities) ,  while  in  the  static  case  only  the  dislocation 


densities  appear  in  the  formulation.  Thus,  higher  order  shape  functions  for 


Au^  are  desirable  in  the  dynamic  case.  The  complementary  conservation 


integral  J,  gives  rise  to  more  singular  BIE's  which  offer  no  advantages  over 


the  conventional  equations. 


The  representation  integral  for  the  traction  components,  Eq.(3.23),  can 


be  used  to  derive  BIE's  for  general  boundary  value  problems  (not  necessary 


5555 


cracks)  of  time -harmonic  elastodynamics  or  elastostatics .  The  advantages 
and  drawbacks  of  this  approach  compared  to  the  conventional  formulation  have 
been  discussed  in  the  paper  by  Hu  (1987) . 

As  pointed  out  by  Hu  (1987),  new  BIE's  can  be  derived  from  other 
conservation  integrals.  Following  essentially  the  same  procedure  as 
described  in  section  3,  the  present  authors  have  obtained  another  set  of 
representation  formulas  for  combinations  of  Uj^(x)  and  u^  j  from  the  well 

known  M  and  integrals  (see  Knowles  and  Sternberg  (1972),  Budiansky  and 

Rice  (1973)).  The  significance  of  these  representation  formulas  and  their 
associated  BIE's  for  solving  boundary  value  problems  is  under  investigation. 
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Appendix  A;  Green's  functions 


The  Green's  function  for  the  3-D  time -harmonic  elastodynamic  state  is 
given  by  (see  Tan,  1975,  Achenbach  et  al . ,  1982) 


exp(ik,jr)  exp(ik^r) 


exp(ik,^r) 

ik  47r/ir  ^ik 


(A.l 


where 


r  - 


(A. 2; 


The  function  u^^(x-y)  denotes  the  displacement • in  the  i-direction  observed 
at  position  x  due  a  unit  force  in  the  k-direction,  applied  at  position 
The  corresponding  components  of  the  stress  tensor  follow  from  Hooke's  law 


C 


G 
u  , 

ijmn  mk,n 


(A. 3 


Similarly,  the  Green's  functions  for  the  2-D  plane  strain  and  anti-plane 
strain  time -harmonic  elastodynamic  states  are 


u°  ( X  -  V ) 
a-y  - 


u^3(x-y) 


[H^^^k,^r)  - 


i—  (k  r) 


H^^^(k-r)]  +  k^S 

o  L  ,0^7  T  a7  o 


H^^^k,^r)j- 


(A. 4; 


(A. 5; 


respectively,  where  denotes  the  Hankel  function  of  the  first  kind 

and  zeroth  order.  Expressions  for  the  stress  components  can  be  obtained  by 
using 


C  uf  ,  for  plane  strain 


q/37 


(A. 6) 


'3a3  -  ^  ^33, a 


for  anti-plane  strain 


(A. 7) 


The  Kelvin  solution  for  the  3-D  elastostatics  case  may  be  written  as 


77 — ^  [(3-4i/)i..  +  r  .r  ,  ] 
l(>irn(l-u)v  '•  ik  ,1  ,k 


while  for  2-D  elastostatics  we  have 


(A. 8) 


u^  (x-y) 

07  -  '‘■ 


7 — yz - -  [(3-Uu)8  inr 

87r/Li(l-i/)  07 


r  r  J 

,a 


(A. 9) 


(A.  10 


for  plane  strain  and  anti-plane  strain,  respectively.  The  corresponding 
stress  components  for  elastostatics  follow  from  the  Eqs.(A.3),  (A. 6)  and 
(A.  7)  . 

We  note  that  both  the  dynamic  and  the  static  Green's  functions  possess 


the  same  singularities  at  r  -  0,  namely, 


“"ijk  ~ 

G  __  1_ 


as  r  -►  0 


(A. 11 


for  the  3-D  case,  and 


1 

u|3 

•v* 

^3q3 

i 

,  5 . 

^"303,^ 

as  r  -♦  0 


(A.  12 


for  2-D  plane  strain  and  anti-plane  strain.  All  derivatives  in  the  Green's 
functions  are  understood  to  be  with  respect  to 

Appendix  B :  Asymptotic  Expansions  of  the 
Elastodvnamic  Green's  Functions 
For  1^1  »  |x|,  the  following  approximation  holds 


r  -  (x-^l  =»  1^1  -  2»x 


(B.l) 


where  ^  denotes  the  unit  vector  along  By  using  (B.l),  asymptotic 
expressions  for  the  3-D  elastodynamic  Green's  functions  are  obtained  as 


g  f  ^  exp(ik^lyl) 

^  exp(-ik^jr.x)  , 


(B.2) 


Fig.  I: 


Curved  crack  of  arbitrary  shape; 


(a)  x^x^-plane. 


(b)  Cop  view. 
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path  independent  integral 
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An  elastodynamic  conservation  integral,  the  integral,  is  employed  to 
derive  boundary  integral  equations  for  crack  configurations,  in  a  direct  and 
natural  wav.  These  equations  immediately  have  lower  order  singularities  than 
the  ones  obtained  in  the  conventional  manner  by  the  use  of  the  Betti-Rayleigh 
reciprocity  relation.  This  is  an  important  advantage  for  the  development  of 
numerical  procedures  for  solving  the  BIE's,  and  for  an  accurate  calculation  of 
the  strains  and  stresses  at  internal  points  close  to  the  crack  faces.  For 
curved  cracks  of  arbitrary  shape  the  BIE's  presented  here  have  simple  forms, 
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and  they  do  not  require  integration  by  parts,  as  in  the  conventional  formulation. 
For  the  dynamic  case,  the  unknown  quantities  are  the  crack  opening  displacements 
and  their  derivatives  (dislocation  densities) ,  while  for  the  static  case  only  the 
dislocation  densities  appear  in  the  formulation.  For  plane  cracks  the  boundary 
integral  equations  reduce  to  the  ones  obtained  by  other  investigators. 
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